Roadmap on Electronic Structure Codes in the Exascale Era

Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing.Comment: Submitted as a roadmap article to Modelling and Simulation in Materials Science and Engineering; Address any correspondence to Vikram Gavini ([email protected]) and Danny Perez ([email protected]

First-principles Investigation of Small Polarons in Metal Oxides

Ein limitierender Faktor der Leitfähigkeit ist die Wechselwirkung der Ladungsträger mit polaren Phononenmoden; das resultierende Quasiteilchen wird als Polaron bezeichnet. Die Stärke der Elektron-Phonon (el-ph)-Wechselwirkung bestimmt die Stärke der Lokalisierung des Polarons, die z.B. die Charakteristik der Temperaturabhängigkeit der Mobliltät definiert. Wir fokussieren uns auf Metalloxide mit starker (el-ph)-Kopplung, bei der sich kleine Polaronen bilden. Die Dichtefunktionaltheorie wird häufig für zur Simulation von Polaronen verwendet. Jedoch treten hierbei zwei Schwierigkeiten auf: Die Sensitivität der berechneten Eigenschaften in Abhängigkeit der Fehler im Austausch-Korrelations (XC)-Funktional und der Effekt der endlichen Superzellgröße. Beide Probleme werden in dieser Arbeit untersucht. Die Polaroneneigenschaften werden auf einer modifizierten Potentialoberfläche (PES) berechnet. Durch Variierung des Anteils der exakten Austauschenergie im hybriden HSE-Funktional zeigen wir, dass das modifizierte PES-Modell deutlich die Abhängigkeit der Polaroneneigenschaften vom XC-Funktional reduziert. Basierend auf dem Potential der el-ph-Kopplung von Pekar leiten wir das korrekte elastische langreichweitige Verhalten des Polarons und darauf aufbauend eine Korrektur für den Fehler durch die endliche Superzellgröße her. Diese Erkenntnisse werden durch ausgiebige Tests an MgO und Rutil TiO2 überprüft. Die oben beschriebene Methode wird zur Untersuchung des Einflusses der Kristallstruktur auf die Bildung von Polaronen in Rutil und Anatas TiO2 und in der β- und κ-Phase von Ga2O3 angewendet. Während in Rutil nur kleine Elektronpolaronen stabil sind, finden wir in Anatas nur stabile Lochpolaronen. Hingegen existieren in beiden Phasen von Ga2O3 nur stabile Lochpolaronen, jedoch mit deutlich unterschiedlichen Bindungsenergien. Dadurch kann durch Verwendung unterschiedlicher Kristallstrukturen Eigenschaften wie Leitfähigkeit und Mobilität der Ladungsträger beeinflusst werden.An important factor limiting the conductivity is the interaction of the charge carrier with polar phonon modes. Such a phonon-dressed charge carrier is called polaron. The strength of the electron-phonon (el-ph) interaction determines the localization of the polaron, which in turn e.g. defines its characteristic temperature dependence for the charge-carrier mobility. We focus on metal oxides with strong el-ph coupling, where small polarons are formed. Density-functional theory is often used for calculating properties of polarons. However, there are two challenges: sensitivity of the calculated properties to the errors in exchange-correlation (XC) treatment and finite-size effects in supercell calculations. In this work, we develop an approach that addresses both challenges. The polaron properties are obtained using a modified neutral potential-energy surface (PES). By changing the fraction of exact exchange in the hybrid HSE functional we show that the modified PES model significantly reduces the dependence of the polaron properties on the XC functional. Based on Pekar's potential for the long-range el-ph coupling, we derive the proper elastic long-range behavior of the polaron and a finite-size correction for the polaron properties. These findings are proofed by an extensively test for rock salt MgO and rutile TiO2. Finally, the approach is used to investigate the influence of the crystal structure on the polaron properties for rutile and anatase TiO2, as well as for the monoclinic β- and orthorhombic κ-phase of Ga2O3. While in rutile TiO2 only small electron polarons are stable, only small hole polarons are found in anatase. Further, small hole polarons exist in both Ga2O3 polymorphs but have significantly different binding energies. Thus, we conclude that growing crystals of the same material but with different structure can be used to manipulate conductivity and charge-carrier mobility

First-Principles Investigation of Small Polarons in Metal Oxides

Ein limitierender Faktor der Leitfähigkeit ist die Wechselwirkung der Ladungsträger mit polaren Phononenmoden; das resultierende Quasiteilchen wird als Polaron bezeichnet. Die Stärke der Elektron-Phonon (el-ph)-Wechselwirkung bestimmt die Stärke der Lokalisierung des Polarons, die z.B. die Charakteristik der Temperaturabhängigkeit der Mobliltät definiert. Wir fokussieren uns auf Metalloxide mit starker (el-ph)-Kopplung, bei der sich kleine Polaronen bilden. Die Dichtefunktionaltheorie wird häufig für zur Simulation von Polaronen verwendet. Jedoch treten hierbei zwei Schwierigkeiten auf: Die Sensitivität der berechneten Eigenschaften in Abhängigkeit der Fehler im Austausch-Korrelations (XC)-Funktional und der Effekt der endlichen Superzellgröße. Beide Probleme werden in dieser Arbeit untersucht. Die Polaroneneigenschaften werden auf einer modifizierten Potentialoberfläche (PES) berechnet. Durch Variierung des Anteils der exakten Austauschenergie im hybriden HSE-Funktional zeigen wir, dass das modifizierte PES-Modell deutlich die Abhängigkeit der Polaroneneigenschaften vom XC-Funktional reduziert. Basierend auf dem Potential der el-ph-Kopplung von Pekar leiten wir das korrekte elastische langreichweitige Verhalten des Polarons und darauf aufbauend eine Korrektur für den Fehler durch die endliche Superzellgröße her. Diese Erkenntnisse werden durch ausgiebige Tests an MgO und Rutil TiO2 überprüft. Die oben beschriebene Methode wird zur Untersuchung des Einflusses der Kristallstruktur auf die Bildung von Polaronen in Rutil und Anatas TiO2 und in der β- und κ-Phase von Ga2O3 angewendet. Während in Rutil nur kleine Elektronpolaronen stabil sind, finden wir in Anatas nur stabile Lochpolaronen. Hingegen existieren in beiden Phasen von Ga2O3 nur stabile Lochpolaronen, jedoch mit deutlich unterschiedlichen Bindungsenergien. Dadurch kann durch Verwendung unterschiedlicher Kristallstrukturen Eigenschaften wie Leitfähigkeit und Mobilität der Ladungsträger beeinflusst werden

First-principles supercell calculations of small polarons with proper account for long-range polarization effects

We present a density functional theory (DFT) based supercell approach for modeling small polarons with proper account for the long-range elastic response of the material. Our analysis of the supercell dependence of the polaron properties (e.g., atomic structure, binding energy, and the polaron level) reveals long-range electrostatic effects and the electron-phonon (el-ph) interaction as the two main contributors. We develop a correction scheme for DFT polaron calculations that significantly reduces the dependence of polaron properties on the DFT exchange-correlation functional and the size of the supercell in the limit of strong el-ph coupling. Using our correction approach, we present accurate all-electron full-potential DFT results for small polarons in rocksalt MgO and rutile TiO2.Peer reviewe

Measurement of the charge asymmetry in top quark pair production in pp collisions at √s = 7 using the ATLAS detector

A measurement of the top-antitop production charge asymmetry A is presented using data corresponding to an integrated luminosity of 1.04 fb −1 of pp collisions at TeV collected by the ATLAS detector at the LHC. Events are selected with a single lepton (electron or muon), missing transverse momentum and at least four jets of which at least one jet is identified as coming from a b -quark. A kinematic fit is used to reconstruct the event topology. After background subtraction, a Bayesian unfolding procedure is performed to correct for acceptance and detector effects. The measured value of A is , consistent with the prediction from the Monte Carlo generator of A =0.006±0.002. Measurements of A in two ranges of invariant mass of the top-antitop pair are also shown